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m THE GTDS FOR 

SAT ELLIT E -to -SAT E ELITE TRACKING DATA 


I. INTRODUCTION 


NASA is currently examining techniques for utilizing a Tracking Data 
Relay Satellite System (TDRSS) to track low orbiting satellites. Such 
a system, employing a relay satellite in a geostationary orbit, has the 
obvious advantage of being able to track satellites over a large portion of 
the globe with only a minimum number of ground stations. However the 
TDRSS concept gives rise to other difficulties not otherwise considered 
in the tracking of satellites from the ground. One such difficulty is in 
obtaining refraction corrections through large portions of the atmosphere. 

In the TDRSS mode the ray path between target and relay satellites 
frequently traverses large distances, and the inhomogeneous nature of the 
ionosphere becomes much more evident, presenting a far greater problem 
raider these circumstances than in ground tracking. One must contend 
with geographic as well as diurnal ionospheric effects, for the line of sight 
between satellites can cross a day -night interface or lie within the equatorial 
ionosphere. 

The purpose of this effort was to examine these various effects and to devise 
a method of computing ionospheric refraction corrections to range and range - 
rate measurements with sufficient accuracy to be of use in orbit determinations. 
The ionospheric model for satellite-to-satellite (SST) refraction corrections 
obviously has to be global in scope. The Bent Ionospheric Model {refs, la, lb, lc) 
is ideally suited for this purpose, particularly as it is already incorporated 
in the Goddard Trajectory Determination System (GTDS) for correcting ground 
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trackmg data. The horizontal electron density gradients encountered along 
the SST ray paths can thus be generated using this model. The vertical 
density profiles predicted by the Bent model are of the form shown in figure (1). 

Making use of the Bent hrr^spher ic Model a method of computing refraction 
corrections through lar o spheric gradients was devised and implemented 

into the GTDS. Herewith report on the various considerations taken in 
designing and implementing this SST refraction correction algorithm. The 
report is structured as follows. Section II is devoted to a discussion of 
satellite to satellite ray path geometries and of their relation to the ionospheric 
model. Sections HI and V describe the development and eval*. 'on of the 
technique used to integrate the horizontal electron density grat .ents encountered 
along a ray path through the ionosphere. Section IV describes the 
implementation of the technique into the GTDS, including program flow 
diagrams contained in Appendices B to F . 


II. SAT ELUTE -T O -SAT E LLIT E RAY PATH GEOMETRIES 


One of the first considerations was to define all the possible ray path 
geometries which could be encountered in the SST mode. This was 
necessary to ensure that all possible ways of traversing the electron density 
model could be accommodated by the algorithm. Since numeric interpolation 
and integration schemes are employed, discontinuities of any kind have to be 
avoided, and pai'ticular care must be exercised whenever the ray paths cross 
domains where the electron density profiles take on different functional 
forms . 

Figure 2 represents three of the general categories of orbits considered: 

(1) target satellite above the ionosphere, (2) target satellite high in the 
ionosphere (above the maximum density layer), and (3) target satellite low 
in the ion isphere (near or below the maximum density layer). Two additional 
tracking configurations are shown, ground tracking of the relay satellite and 
the target satellite. The ground tracking case must also be included, because 
the ground to relay satellite leg is an integral part of the SST tracking 
geometry. The important features to note about these various geometries 
are discussed below. In addition, five specific satellite to satellite orientations 
are taken up in greater detail, including program flow diagrams, in section 
IY and Appendices B-F. Figure 3 shows some of the parameters, to be 
encountered in the discussion below, which play a significant role in defining 
the electron density gradient function and in the integration of the densities along 
the ray path. Definitions of other terms and parameters used in this report 
can be found in the glossary (Appendix H). 

A. Target Satellite Above the Ionosphere 

Tills class of orbits gives rise to the longest transionospheric ray paths to be 
encountered. The geometry involved is shown in figure 2 and figure 3 . 

Moving from the target to the relay satellite the ionosphere is entered at point 


-4- 


P1 and left at point 1 Midway between PI and P2 lies PC, the point of 
closest approach to the center of the earth. Tracing the target satellite in 
its orbit from the time it completes occupation, the following changes in 
the segment (P1,P2) are observed. At first, the ray path enters the 
ionosphere at PI, encounters a density peak at HM, exits through the bottom 
of the ionosphere at PB, reenters the ionosphere at a second point PB, 
encounters another density peak at HM and finally exits once more at P2. 

In effect two distinct ionospheres are encountered along the ray path. It is 
important that the density peaks at HM along the ray path be properly placed 
and the density correctly evaluated. Errors made at these locations can be 
significant, for the densities here are orders of magnitude larger than 
at the top and bottom of the ionosphere. 

As the target moves on in its orbit, PC rises higher above the surface until . 
it crosses HMYM, the bottom of the ionosphere. At this time, the segment 
(PI, P2) lies entirely in the ionosphere, but two density peaks are still 
encountered. Care must be exercised to compute the:.-- maximum densities 
properly in order to minimize possible errors in integrating the total 
electron content along the path. 

Next, PC approaches the height of maximum density HM. In this case 
a significant portion of (P1,P2) lies within the peak density region and special 
attention must be paid not only to obtaining the maximum density, but also to 
calculating accurate density gradients. 

As PC rises further, segment (P1,P2) becomes shorter, and at the same time 
the electron density decreases exponentially. The net effect is a very rapid 
decrease in integrated electron content until the satellite to satellite ray path 
leaves the ionosphere entirely. It is difficult to generalize the behavior of 
the range and range rate corrections for satellite to satellite geometries 
because of the many variable ionospheric conditions that can be encountered 
(day to night gradients, latitudinal variations , equatorial anomaly, variable 
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h eight cf maximum density, etc, ), However, Tor those cases when the target 
satellite is above the ionosphere, the corrections do follow a certain pattern. 

As the target merges from occultation, the range correction increases from 
a non-zero value to a maximum, achieved when PC crosses the peak density 
at HM. From this maximum value, the range correction then decreases 
monotonically to zero. The range rate is even more difficult to generalize, 
for it is obtained from the time rate of change of the range correction. It 
attains a maximum value before PC reaches HM, a minimum value as PC rises 
into the exponential layers of the ionosphere, and goes to zero as PC approaches 
the top of the ionosphere. The sequence is repeated, (in reverse with 
sign changes on the range rate corrections), as the target satellite approaches 
occultation on the other side of its pass. 

B. Target Satellite High In The Ionosphere ( above peak density) 

This geometry is similar to the previous case of having a satellite above the 
ionosphere, with the difference that the range correction never goes to zero. 
When the target satellite emerges from occultation, the range correction 
increases until PC crosses HM. The correction then declines to a minimum, 
which is achieved when the target passes beneath the relay satellite. As the 
target continues in its orbit a second range maximum is reached when PC 
again crosses HM. The magnitude of the peaks will depend on factors such 
as the height and inclination of the orbit of the target and, of course, on the 
ionospheric state. If the range correction peaks occur on opposite sides of 
a day-night interface then they could differ in magnitude by a factor of ten. 
Under such circumstances it is not possible to discuss the general appearance 
of the range rate corrections. 

C. Target Satellite Low In The Ionosphere - (below peak density) 

The SST configurations, in which the target is in an orbit below the peak 
density of the ionosphere, result in the most complex range and range rate 
correction signatures. The refraction correction is determined by two factors. 
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One factor is the ge metric length of the SST ray path through the ionosphere, 
and the other is the ionospheric state. The geometric effect is such that 
maximum path lengths are achieved near occulta tion when the point of closest 
approach, PC, is still above the bottom of the ionosphere, HMYM, (fig. 3). 

In the presence of a homogeneous ionosphere, these positions coincide with 
the position of the range correction maxima. The maxima would be 
■v, symmetrically situated about the midpoint of each pass of the target. The real 

ionosphere, however, is inhomogeneous and the large scale gradients in 
this complex medium will dominate the geometric effects of increased path 
length. Hence the presence of a day to night interface or the crossing of the 
equatorial ionosphere (in the course of a satellite pass ), will determine 
the magnitude and placement of the range correction maxima (or maximum). 
Typical ionospheric gradients are illustrated in figures 4. For equatorial target 
satellite orbits, the day-night (sunrise or sunset) gradients will be of 
greatest significance, while for polar orbits the crossing of the equatorial 
ionosphere will be the dominant ionospheric effect. Orbits of intermediate 
inclination will be affected by a combination of both of these large gradients 
and the resultant range corrections will have the most complex signatures of all 
SST configurations. The range rate corrections are impossible to 
discuss in general terms for these cases. 

D. Direct Tracking of Satellites From the Ground 

This ray path configuration completes the satellite to satellite tracking loop 
by linking the relay satellite to the ground station. The SST refraction 
algorithm developed for the GTDS therefore has the additional capability of 
computing refraction corrections for direct tracking of a target satellite 
from the ground. When tracking the relay satellite, the ionospheric range 
correction to the relay will be nearly constant during the course of one pass 
of the target. Since the relay satellite is stationary relative to the ground 
station, the only variation in the range correction will be due to diurnal 
variations of the ionosphere. 


For direct tracking of the target, the situation is quite different, for the motion 
of the satellite adds a significant factor to the range rate correction. Combining 
geometric and ionospheric effects one can obtain a variation in the range 
correction of a factor of 30 just in the visible sky over one tracking station 
(ref, 6). The variation in range correction will depend on the orbit of the 
target as well as on the ionospheric state, just as it does for SST configurations. 
The dominant factor in the range rate correction under these conditions is 
the velocity of the Carget relative to the ground station. In general the rate 
correction should achieve a positive maximum as the satellite rises in the 
sky over the station, and a negative maximum as the satellite sets over the 
station. The relative magnitude and placement of these maxima depends on 
the ionospheric state. 


III. INTEGRATING TOTAL. ELECTRON CONTENT ALONG SATELLITE 
TO SATELLITE RAY PATHS 


The major problem in computing range and range rate corrections lies in 
accurately defining the horizontal electron density gradients which are 
encountered along the large distances characteristic of satellite to satellite 
(SST) geometries. A satellite at 1000 km can achieve SST ray path which 
will subtend earth central angles of 80° to 90° totally immersed in the 
ionosphere. Figure 4a gives some idea of the magnitude of the gradients. 

It is a world map of the ionospheric critical frequency fOF2 projected onto 
a shpere at height 300 km above the surface of the earth. fOF2 is related 
to the maximum electron density by the equation 

N, =1.24x 10 10 x (fOF2) 2 el/m 3 (1) 

where fOF2 is in MHz. Figure 4b is a world map of the height of the maximum 
electron density HM. It too shows a highly complicated structure. Both of 
these parameters play a critical role in predicting the vertical density profile 
using the Bent model. In addition to providing N B , f 0F2 is also one of the 
key parameters used to compute the thickness of the parabolic regions of the 
profile, and the decay constants in the exponential layers (figure 1 ). 

Ir.-iig rating the profile vertically gives the vertical electron content 

btop 

N v = f N (h) dh (2) 

h «-y. 


Figure 4c shows how this quantity varies globally on a typical day. The 
extent to which the ionosphere is inhomogeneous becomes strikingly obvious 
in this figure. 


A. The Gradient Function 


The total electron content between any two points , PI and P2, in the 
ionosphere is obtained by integrating 

N t = 

PI 

where N(s) is the electron density function, giving the electron density at 
position s- s(x, y, ss) or s( X, 0 , h) along the path.and 0,h are the longitude, 
latitude, and height above the surface, respectively. Since N(s) does not 
have a simple analytic representation some new means had to be devised to 
perform the integration in eq. 3. Accordingly we chose a quadrature 
integration teclmique. Using quadrature integration it is not necessary 
to know the exact functional form of the integrand. It suffices to have values 
of the integrand at specific "quadrature" points. 

The problem of dealing with a complicated function like N(s) was then reduced 

to computing N{s) at some specific points along the path. To simplify these 

computations N(s) was separated into 2 parts. The approach was to consider 

the ionosphere as a three dimensional medium having a vertical electron 

density distribution given by the profile in figure 1. If the ionosphere were 

homogeneous then one such vertical density profile would suffice to define the 

ionosphere any place on the globe. Iicwever, the ionosphere is not homogeneous 

therefore, each location on the globe has its own profile. But, since the profiles 

all have the same basic shape it is possible to separate this basic height 

dependence from the positional variations. The basic height dependence is 

represented by a reference profile N re f { ^ j. , §. e £ , h ). The positional 

dependence is then reduc ed to perturbations or deviations from the reference 

profile which are represented mathematically by a function G(s), called 

the gradient function. ALong any ray path the line of sight electron density 

distribution N(s) is then given by one reference profile modified by a 
gradient function, as in eq.4. 


| N (s) ds (3) 


Figure 5 may be helpful in illustrating the relationship between the functions 
N,N re £ and G. It now remains to evaluate the gradient function G which, 
b equation 4, is 


G(s) 


JM 


^ref 


U ,, e 

ref ref 
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B. The Deviation Function 


During the course of numerical evaluations of the function G it became expedient 
to represent the gradient function as 

G(s) = 1 + 6 (s) (6) 


In this form 6(s) isolates the deviation of the density at s from the 
reference profile density. For this reason 6 (s) is called the deviation 
function. Now, by equations 5 and 6 


6<s) 


_ N (s) 

^ef <kf ' ref * h > 


(7a) 


Equation 7a is the defining equation of the deviation function. 


If the variable s is replaced by spherical positional parameters 8, and 
h, it may become clearer how 6 is evaluated in a practical application. 


6 (s) 


n (\,e,h) 


N i\ , 0 , h 

re l rel ref 


(7b) 
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fe (s) is defined at several discrete points {up to 22 points in the SST 
formulation ) between the ends of the ray path, thereby producing a table 
of values, fe (s^ ), where 


6 . 

1 



N. (-X., 9. ,h.) 

1 l 1 1 

N . (X ,0 , h.) 

lef, i ref ref i' 


-1 


(7c) 


Theoretically any reference profile could be used to evaluate the 6. 

in eq.7c. However, when defining 6 (s) using only a few discrete points, 

it was found necessary (in order to maximize numerical accuracy) to choose 

the reference profile at the point (\ , 0 , IIM) along the ray path where the 

maximum electron density is encountered. All the profile parameters are 

generated at the location (X , , 0 , ) using the Bent Ionospheric Model. 

rex ref 

They are later used, as needed, to compute the densities N ^ ^ at the 
heights h. . Once-the reference profile has been chosen, the densities N. 
are computed. To do this a complete ionospheric profile, i.e. a complete 
set of profile parameters, is predicted by the Bent Model, at each location 
( ^i ’ ®i ) in turn# Using the profile parameters, the density hL at each 
point (X. i , 0. , In ) is obtained by substituting the height In into the equations 
definmg the profile. The profile equations are listed on figure 1 . 

To obtain values of 6 at intermediate points of the 6,{s.) table, several 

ix 

numerical techniques were investigated. The first approach was to represent 
fe by a simple polynomial. Based on examinations of world maps of fOF2 
and N (figures 4a and 4c) it appeared that a polynomial of 1st or 2nd 
order would suit the mid and high latitude regions of the ionosphere while a 
3rd or 4th order polynomial would fit the equatorial regions. Accordingly, 
fewas computed at numerous points along various ray paths and polynomials 
were fitted to these data in a least squares sense. It was necessary to go to 
sixth order fits before numerical results became stable enough to perform 
range differencing. Fils of this order would require too many data points, 
hence, the polynomial representation of 6 was 


abandonned, and instead, numerical interpolation schemes were investigated. 
Interpolation using a cubic spLine technique (ref. 3) was finally chosen as 
the most satisfactory means of defining 6 at intermediate points in the 
6. (s.) table. Details concerning the number of points needed to 

accurately define 6. as well as the placement of those points, are discussed i 
section V below. 

As previously mentioned, the reference profile is always chosen at the 
point of highest density along a ray path. For direct tracking of satellites 
from the ground, this point occurs when the ray path crosses HM. In the SST 
mode it is possible for the ray path to cross HM twice. In such cases, the 
ray path is divided into two legs, the division being made at i-the point of 
closest approach, PC . A separate reference profile and deviation function 
are then obtained for each leg. For the cases when PC is above HM, the 
reference profile is chosen at PC, for, being the lowest point in the ionospher 
PC becomes the highest density point on the ray path. For those geometries 
where the target satellite itself is the lowest point in the ionosphere, the 
refer ence profile is chosen at the target satellite. 

This concludes the description of how the electron density is computed, at any 
point along an integration path through the ionosphere, as the product of a 
reference profile and a gradient function. Both the reference profile and 
the gradient function have been defined and discussed. Let us now turn to the 
problem of actually integrating the electron density N(s) from point Pi to 
point P2, as required by equation 3. 

C . Defining The Integration Path 

The integration of the total electron content now proceeds according to 
equation 8, obtained by substituting eqs. 6 and 4 into eq. 3. 

P2 

N t=J (N ref * (1+5) ) ds (8) 

PI 




The integration in eq.8 is executed assuming the path (Pl,P2) ie a straight line 
(refer to figure 3 for definition of points). To assume straight line radio 
propagation means that e.ffects of refractive bending are being neglected. 

This assumption is valid for satellite to satellite tracking, since at S -band 
frequencies refractive bending is negligible. 

The integi ion interval (P1,P2) is not always treated as a single interval. 

The interval is sometimes split into two intervals (PC, PI) and (PC,P2) 
and the integrations then summed. The distinction between these two modes 
of integration is made on the basis of the position of the point of closest 
approach, PC, relative to the height of maximum density HM. If PC 
is above HM, then only one deviation function is defined for the entire path 
(with reference profile at PC) and the interval (PI, P2) is not split at PC. 

If PC falls below HM. and PC lies between PI and P2, the path is divided 
at PC into the two intervals (PC, PI) and (PC,P2). Each of these two intervals 
has its own deviation function, defined with respect to a reference profile 
computed at the intersection of the ray path with HM. In the event that PC 
falls below the ionosphere, the intervals become (PB,P1) and (PB,F2). 



For those SST configurations in which the target satellite lies on the near ; 

side of the earth, as seen from the relay satellite, i. e. PC lies outside 
the interval (PI, P2) (see figure 3), then the interval (PI, P2) is not divided. 

I 

A single deviation function is computed in these cases, relative to a reference i 

profile chosen either at HM, or at the target, if the target is above HM. i 


The purpose in distinguishing between these various ray path configurations, 
is to ensure an accurate definition of the density gradients along the path 
(F1,F2), particularly when the path traverses the region around the peak 
electron density at HM. 


I n performing an actual integration, the interval (P1,P2) is further subdivided 
into segments, so chosen as to reflect the boundaries of the seven layers of the 
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reference profile (fig. 1). The end points of these segments correspond to the 
layer heights HMYM (bottom of the ionosphere}, HM, Hu, Hi, Kg, 1000 km, 
2000 km, and 3500 km. This segmentation of the integration interval is 
necessitated by the requirement of quadrature integration that the function 
to be integrated and its derivatives be continuous over the integration interval. 
To meet these requirements, the integration of equation 8 by quadrature must 
therefore proceed in steps through the layers over which the reference profile 
does not change analytical form. 


D. Numerical Integration of the Total Electron Content 


The Gaussian quadrature method (ref . 4) was selected to evaluate the total 
electron content according to equation 8. Gaussian integration is performed 
as follows: 



b - a 
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( 10 ) 


= Zeros of the Legendre Polynomials (x) 


weights w. - Jf^JT [*>' W] 


_ a J 2ll + 1 fa | J 

remainder R n = ( 2 r,+l/ O 2 ") ' J 
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Tables of the values of x. and w. can be found in many standard texts on 
numerical techniques, including ref. 4 . 

Equation 8 can now be rewritten according to equation 9 to produce the form 
which is implemented in the GTDS for SST refraction correction calculations. 
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N . (s) (1+6 (s) ) ds 
ret 



( 12 ) 


and S = position along the ray path (Pa,Pb) which intersects the height 

j 

of the bottom of the j th layer of the reference profile 

A. = index of the lowest profile layer encountered along the interval 
(Pa, Pb) 

m = index of the highest profile layer encountered along the interval 
(Pa, Pb) 

= position corresponding to lowest height attained by path (Pa, Pb) 
in traversing the ionosphere 

S = position corresponding to highest height attained by path (Pa,Pb) 
m 

in traversing the ionosphere 

Pa,Pb = limits of integration corresponding to integration intervals 

(PI, P2), (PC,P2), (PC, PI), (PB, PI), (PB, P2), as appropriate 
for each SST ray path 

N = reference profile 
ref 

6 = deviation function defining density gradients between points 
Pa and Pb in the ionosphere 


n = number of quadrature points used to integrate electron density 
in j*k layer; n- = 4 for j = l and n:=3 for j=2, 3, ... 7 

Equation 11 is the form in which integrations are carried out for SST 
configurations. The first step in evaluating equation 11 is to determine 
the limits of the j - sum, i.e. the summation over layers, as set by the 
configuration geometry. This procedure was described in section III C, above. 
I ntegration then proceeds from layer l to layer m. 

The quadrature weights and points, w^ and x^, appropriate for each layer are 

obtained from tables. The minimum number of points, m, required to 

accurately integrate N (1+6 ) through each layer, was determined in a 

r cl 

series of evaluation runs in which quadrature integrations were compared 
with ray tracing calculations. In the course of these evaluations it was 
found that at least 4 quadrature points were needed to integrate through the 
first (bi-parabolic) layer of the ionosphere, and at least 3 quadrature points 
were needed for each of the six topside layers (one parabolic and five 
exponential layers). 


Using the appropriate quadrature points for the layer lying between points 
Sj and along the ray path, equation 12 provides the positions s^ at which 

the deviation function 6 and the reference electron density N. are to be 
evaluated. The reference density is obtained by substituting the height h| 
corresponding to the point s^ (i.e. s^=s ( 0^, h. ) ) into the ionospheric 

profile equations (fig. 1) containing the reference profile parameters. The 
deviation function is interpolated at s. using a cubic spline interpolation as 
described in section III-B, above. The products w^ * N (sj ) * (1+6 (s|) ) 
are then formed and summed, and the totals for each layer added according 
to equation 11 . 


If the integration interval (PI, P2) (see figure 2) has been split into two 
intervals, as described in section IH~C, then this summation procedure is 
repeated seperately for each interval and the results are then added to produce 


a total electron content, N T , between points PI and P2. This is the basic 
quantity which is used to compute the ionospheric refraction correction to 
range measurements. As described in section IV, a simple scaling of N| 
by the tracking frequency produces the range correction. 

Details of some of the evaluations undertaken during the development of this 
refraction correction technique are discussed in section V below. These 
include studies of the definition of the deviation function, the selection of 
the integration paths, and the definition of the minimum number of points 
required to perform accurate quadrature integration of equation 8. 


IV. RANGE AND RANGE RATE CORRECTIONS 


A. Derivation of the Ionospheric Refraction Corrections to Range and 
Range Rate 


The refraction corrections to range and range rate are obtained according 
to the following equations. 


„ 40. 3 

Range: AR =— N T 

Range-rate: AR = A , t ) . : .. A R W 

A t 


(13) 


(14) 


where the total content N T is defined in equation (3) and evaluated according to 
equation 8, f is the tracking frequency in Hz, At is in seconds and AR and 
AR are in meters and meters /sec . , respectively. The range and the range 
rate corrections so defined have opposite effects on their respective 
observable quantities because of the different manner of propagation of the 
signals from which the ranges and rates are extracted. 

Ranges are obtained from the modulation of the carrier wave of the tracking 
signal, whereas range rates are deduced from doppler counts (phase measure- 
ments). The modulated carrier propagates as a wave group through the 

ionosphere and its behavior is characterized by the group index of refraction n 

8 

Since the ionosphere is a dispersive medium the group and phase (n^) indices 
of refraction are not equal. Hence, phase propagation must be characterized 
by the phase index of refraction. Fortunately, at S-band frequencies n„ and 

O 

Up are simply related (ref. 2) by 


For satellite tracking frequencies, the Appleton- Hartree equation (ref. 2) 
simplifies to 


"p = 1 ~ k “P~ OM 

where k= (e® /4 tt 2 tf 0 m) =80. 5 and f is in Hz. By using the binomial 
expansion and equation (15), we get 

n g = 1+k 73“ U6b) 

Using the binomial expansion again and substituting for k, eqs. (16) 
reduce to the more familiar 


40. 3 N 

(17a) 

“p = 1 - 

. , 40.3N 


n„ = 1 + — 

g f 2 

(17b) 


One can now define the radio path R of a signal propagating through the 
ionosphere (figure 6) 



ds . 


(18) 


Substituting eqs. (17) 
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R 


P 2 „ _ 40. 3 N v . 

= (! + — J5 ) ds 


= R, 


- 40.3 
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N ds 


(19) 
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where R 0 is the straight Line range in the absence of any ionosphere. Equation ( 19 ) 
defines the range correction &R (equation 13), such that 

R=R 0 + AR 

Hence for phase and group measurements the radio paths are 


R u 
phase 

= R 0 — AR 

(20a) 

R 

group 

= R 0 + A R 

(20b) 


Range measurements are then corrected according to 


Range , = Ranee + AR 

measured vkcuum 


(21b) 


while range - rates (phase measurements) according to the time derivative 


of eq. (20a) 


Range-rate . = Range-rate — AR 
measured vacuum 


(21a) 


In applying corrections to observations in the GTDS, the convention is such 
that the corrections are always subtracted from the observations, therefore 


Ranee-rate , , =Range-rate . — (~AR)=Range-rate (22a) 

measured, corrected “measured vacuum 


Range , , . =Range , — A R=Range 

measured, corrected measured vacuum 


(22b) 


Equation (22a) shows why the sign on the range rate correction is always 
changed before it is applied to the tracking observation. 


An excellent discussion of this phenomenon from an electrical engineering 
standpoint can be found in ref. lc . 

B. Correcting the SST Data 

Integration of the total electron content encountered along SST ray paths 
through the ionosphere is controlled in the GTDS by subroutine IONSST. 

Th"' selection of integration paths and the definitions of gradient functions is 
all performed in this subroutine. Program flow diagrams of IONSST can be 
found in Appendix B. 

For illustration purposes, five test cases have been traced through the flow 
diagrams of IONSST for various SST geometries. The flow of the logic for 
each case is traced by a heavy line in Appendices B through F. Included are 
three cases involving a single integration path but with differently defined 
reference profiles, one case involving two integration paths requiring two 
separate reference profiles, and one case depicting direct tracking from the 
ground. The geometries are illustrated in figures 7 through 11, and the 
significant points to note for each case are summarized below. 

Case 1 (figure 7) - The target satellite is in the topside ionosphere above HM, 
and is located between PC and P2. The deviation function is defined relative to 
a reference profile at PI, and from 2 to 9 points are used to define the function 
over the interval (P1,P2). The exact number of points depends on the height 
of the point of closest approach, PC, and on the length of the path (P1,P2). 

The closer PC approaches HM, the larger the number of points used to specify 
the deviation function. The integration interval (Pi,P2) is divided into 
segments corresponding to the topside layers of the reference profile, and 
three-point Gaussian quadrature is used integrate the electron content 
through each of these segments. The flow diagram for case 1 is in Appendix B. 




1 


Case 2 (figure 8) - The target satellite is above the ionosphere. PC lies 
above HM so only one reference profile (at PC) is used to define the 
deviation function between PI and P2. Again, depending on the height of PC 
and the length of the path (P1.P2), anywhere from 3 to 17 points, spaced at 
intervals of equal earth central angle, are used to specify the deviation 
function over the interval (P1,P2). Integration then proceeds in layer 
segment from Pi to P2 using tnis deviation function and the reference profile. 
The integration through each segment is executed by a 3-point Gaussian 
quadrature. The flow diagram for case 2 is in Appendix C. 


Case 3 (figure 9) - The target satellite is in the bottomside ionosphere 
below the maximum density layer HM. The position of HM along the path 
(Pi,P2) is obtained by an iterative process. This being the maximum density 
encountered along. the path, the reference profile is computed at this location. 
Next, a profile is computed at PI in order to define the deviation function 
at this point. A comparison is made between the bottomside ionospheres 
at PI and at the reference profile. If it is found that the ionosphere at PI 
extends to a lower height than at the reference point HM, then the bottomside 
of the reference profile is modified to reflect this difference. This is 
achieved by setting the reference profile parameter YM (ref) to the value of 
YM at PI. This precaution is taken to avoid divisions by zero when defining the 
bottom of the ionosphere. The deviation function definition is then completed 
using up to 9 points spaced at equal earth central angles along (P1,P2). An 
additional point is placed halfway between PI and HM to improve numerical 
accuracy. In the event that the earth central angle subtended by (P1,P2) is 
less than three degrees, gradients are neglected and the deviation function is 
set to zero. Integration proceeds from PI to P2 in layer segments. The 
first segment from PI to HM is integrated using 4 quadrature points, while 
the six topside segments are integrated with 3 quadrature points. The flow 
diagram for case 3 is in Appendix D. 


Case 4 (figure 10) - The target satellite is within the ionosphere in a position 

near occupation. PC falls below the bottom of the ionosphere, HMYM, so 

that the ray path exits and reenters the ionosphere at the points PB. These points 
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nii' * be solved for iteratively- The path (Pl,P2) is now subd'vided into 
2 intervals, (PB,P1) and (PB.P2). Reference profiles for each of these 
intervals are computed E- the points HM, which must also be found by iteration. 
Tests are made to determine if the bottomside reference profiles need to be 
extended relative to the bottomside at the points PB. This is to ensure 
properly d* fined deviation functions in the bottomside ionosphere, as discussed 
in the example of case 3. A separate deviation function is defined for each 
interval using up to 9 points, spaced at equal earth central angles. An additional 
point is placed midway between PB and HM. Two integrations are then 
performed, layer by layer, from PB to PI, and from PB to P2, The 
results of the separate integrations are added to produce a net total electron 
content for the path (P1,P2). The flow diagram for case 4 is in Appendix E. 

Case 5 (figure 11) -* Direct tracking of a satellite from the ground. The 
integration interval is (PB.P2). If the satellite is above the ionosphere, P2 
corresponds to the top of the ionosphere at 3500 km. If the satellite is within 
the ionosphere, P2 coincides with the satellite position. The intersections 
at PB and at HM are found iteratively. The reference profile is computed at 
HM along the interval (PB,P2),and from 2 to 7 points are used to define the 
deviation function over this interval. For high elevation angles when the earth 
central angle subtended by (PB,P2) is less than 3° , gradients are neglected 
and the deviation function Ls set to zero. Integration now proceeds by 
quadrature from the first layer up to a maximum of seven layers. The flow 
diagram for case 5 is in Appendix F. 



Once the integration has been completed a return is made to subroutine ) 

1 

CORSST. It is here that the one-way range and range rate corrections are 

j 

computed and applied to the observations, as specified by equations 23,24 and 

2 5. 


(23) 


DMR = 


40. 3 N , (relay) 


( f -V 

up 


_1 
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dwn 


7T ) 
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40. i N r (SST) 

f a 

SST 


3 




= DMR (T+PELT) -DMR(Y) 

BELT 

OM2 = OM2 - DMR/1000 (Range Observations) 

OM2 = OM2 - DRATE/1000 (Range rate Observations) 


( 24 ) 


(25a) 

(25b) 


where N T (relay) = total electron content for the ground to relay 
satellite link 


N t (SST) = Total electron content for the relay to target satellite 
link 


*up * ci ' fesT ~ u P^ n ^» downlink, and satellite to satellite 
tracking frequencies. 

DMR =range correction in meters 
DRATE = range l'ate correction in meters /sec 

OM2 = tracking data observation (GTDS notation) in km of range 
or km /sec of range rate. 

The program flow diagram for subroutine CORSST is found in Appendix A. 

Since it takes two successive range observations to compute a range difference , 
the range rate correction is set to zero if the first observation in any pass 
is a range rate. Furthermore, if the time interval, DELT, between successive 
range rate observations is too large ( 10 minutes), DRATE is again set to 

zero and the range differencing is reinitialized, as for the beginning of a new 
pass. For those cases when the time interval is too small ( < 1 second), 

DRATE is not evaluated, and instead, the last computed value of DRATE 


C. Accuracy of the Corrections 


The ionospheric refraction correction algorithm for satellite to satellite 
tracking data described in this report has not yet been evaluated. Test runs 
using tracking data from the ATS -6 to GEOS -3 SST configuration indicate 
that the algorithm has been successfully integrated into the GTDS. The 
refraction corrections computed by the algorithm, for this configuration, 
range from 0. 5 to 10 meters for range measurements and from 0. 1 to 1. 0 
cm /sec for range rates, at frequencies near 2 GHz. Although they are of the 
right order of magnitude, the corrections appear to be lower than expected 
by perhaps a factor of 2 or 3 in the range correction. Some numeric 
inconsistencies have also been noted, between the integrations in the ground 
to relay link and the relay to target link. Hence, definitive conclusions 
regarding the accuracy of the SST refraction correction algorithm must be put 
off pending further testing and evaluations. 

Some comments can, however, be made regarding the expected accuracy of 
the refraction corrections. The approach of using a reference profile, in 
conjunction with a gradient function, to represent electron density variations 
between two points in the ionosphere, was thoroughly tested in a ground 
tracking configuration for a satellite above the ionosphere. Though this did 
not provide a complete test of the technique, the ground tracking mode did 
allow an in-depth evaluation of ail the requisite numerical techniques including 
quadrature integration through the ionospheric profile, definition of the deviation 
function and spline interpolation of the deviation function. Furthermore, some 
similarities do exist between the two modes. For example, the longest 
path through the ionosphere for ground tracking occurs at 0° elevation when the 
path subtends an earth central angle of about 40° . The longest path in the 
SST mode is about 5 0° , as measured from the point of closest approach to 
the top of the ionosphere. Although the total path in the SST case is 100° , 
for integration purpose it would be treated as two 50° intervals. Hence, the 
low elevation angle lines of sight in the ground tracking mode can provide a 
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fair approximation of SST geometries. 

Range corrections in the ground tracking tests were found to agree with ray 
trace results to within 1% {though with some exceptions) when the deviation 
function was defined by as few as 7 points. The limiting accuracy in 
computing the corrections was entirely dependent on the definition of the 
deviation function, and not on the quadrature integration or spline interpolation. 
Hence, with deviation functions defined by using a few more points (up to 9 or 10), 
it should be possible to achieve the same 1% accuracy in r&nge correction for 
any SST configuration. 

The accuracy of the range rate corrections will vary with the accuracy of the 
range corrections. Experience with ground tracking geometries indicates 
that for accuracies of 1% in the range corrections the range rates should be- 
good to within 2% to 4% . Confirmation of this estimate will have to await 
further testing. 

A more detailed discussion of the ground based evaluation runs can be found 
in section V of this report. 
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V. EVALUATION OF THE SST REFRACTION CORRECTION ALGORITHM 

To date, the refraction correction algorithm as it is implemented in the 
GTDS, has not been evaluated. However, the numerical techniques employed 
in the algorithm were thoroughly tested in the course of program devel- 
opment. The tests were made using a configuration in which a ground sta- 
tion tracks a satellite in orbit above the ionosphere (ref. 7). At low elevation 
angles, ground to satellite ray paths in the ionosphere can be as large 
as 40° of earth central angle. Though this is not as long as some of 
the SST ray paths, these low elevation cases traverse enough ionosphere 
to encounter substantial density gradients. As a result it. was possible to 
confirm the validity of representing the variations in electron density in 
the ionosphere as the product of a reference profile and a gradient function. 

The ground station at which the evaluation tests were made, was located 
at Kashima, Japan (lat; = 35.®95, long = 140?7 ) and the evaluations were run 
for 12h U. T. on 15 January I 969 . The date was chosen to coincide with 
solar maximum, and the location was selected because of its proximity to 
the geomagnetic equator where electron density gradients are high. The 
ionospheric state for this date and time is illustrated in the world maps 
of figures 4. Integrations of total electron content were carried out by 
ray tracing along lines of sight selected at 30 0 intervals in azimuth and for 
19 elevations. From 0° to 10° the elevations were incremented in steps 
of 1°, from 10° to 30° the elevations were incremented in steps of 5 0 , and 
from 30° to 90° in steps of 15“ . 

At low elevations, up to 150 density points were computed to define the elec- 
tron density variations along each line of sight. The points were generated 
at intervals of 20 km along the path, from the bottom of the ionosphere to a 
height of 1000 km, then at 200 km intervals from 1000 km to 3500 km. At this 
height the ionosphere was cut off. As the line of sight approached vertical 
incidence, the total number of generated points decreased to approximately 60. 


The ray tracing program (ref. 5) is a double precision arithmetic application 
of Snell's Law in spherical form. Numerical integration of the equations, 
relating ray bending and electron content to changing local elevation angle 
and refractivity, was performed layer by layer through the ionosphere to 
provide the total electron content. The integration steps were 1 km in 
height below 1000 km and 5 km for heights above 1000 km. Integration of 
the density through each step was carried out using a 3 -point Simpson's 
rule. Electron densities at positions in between the 60 to 150 precomputed 
density points were provided to the ray trace program by interpolation. 
Below the maximum density the interpolation was of fourth order. Above 
maximum density the interpolation was quadratic in the parabolic layer and 
logarithmic in the exponential layers. Comparisons with analytic integration 
at vertical incidence indicated the ray trace total electron content was 

Q 

accurate to 1 part in 10°. 

The ray trace results served as benchmarks for evaluating the gradient 
function technique. The main items of interest were to determine the 
minimum number of points needed for accurate Gaussian integration through 
the gradients and to optimize the calculation of the deviation function. 

A. Optimizing the Tabulation of the Deviation Function 


The deviation function 6 defin ing the electron density gradients between 
two points PI and P2 in the ionosphere is computed at equally spaced 
intervals of earth central angle between the two points. The number of 
intervals to be used is obtained using the empirically determined function, 
NOPF. 


NOPF = 


48x1 0 5 
H+ 18x10 s 


+ 


ECA 
7. 5 


where ECA is the earth central angle (in degrees) subtended by the interval 
(P1,P2), and H ( in meters) is the height of closest approach to the surface of 
the earth of the line segment (P1,P2). The function is integerized in the 
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computer so that it returns a fractionless integer. The definition of NOPF 
is such that more weight is placed on the lower regions of the ionosphere, 
where the gradients are expected to be the most significant. Hence, for 
ray paths crossing through or near the peak density, NOPF returns its 
maximum value of 8 intervals. For the cases when the point of closest 
approach, PC, lies between PI and P2 (see figure 3), the ECA is measured 
from PC, and NOPF is then separately determined for (PC,Pl ) and (PC,P2). 

This means that as many as 16 intervals may be used in defining 6 along 
some of the longer transionospheric ray paths. Up to two more points 
may be added to these, as follows: (1) the point HM, if it lies on the ray path, 
and (2) a point halfway between HM and the lowest point encountered along 
the path. It is therefore possible to have a 6 defined by as many as 22 
points (counting end points). 

The function NOPF was determined on the basis of comparisons with the 
ground based ray tracing calculations. For each line of sight in the benchmark 
cases described above, the deviation function was computed at fifteen 
heights. In each case, the reference profile was at the peak electron density. 
Different combinations of the heights, varying from as few as 6 to all 15 points, 
were then used to define the deviations along the lines of sight. Quadrature 
integration produced a set of range corrections which could be compared to 
the exact ray trace calculations along each ray path. Some of the results of 
these comparisons are presented in Tables 1,2, and 3. The labels in the 
tables mean the following: 

analytical - anaLytic integration thru the ionosphere assuming 
a homogeneous profile at HM 

numerical - Gaussian quadrature integration through the 

ionosphere assuming a homogeneous profile at HM 

numerical - (least squares) - quadrature integration with 

gradients given by 6th order polynomial fitted by | 

least squares j 

numerical - (spline) - quadrature integration with gradients i 

given by cubic spline interpolation through different J 

sets of points j 

■ ? 
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The range and range difference corrections are in meters. Angle corrections 
are in mrad. The percentage difference for each entry is taken relative to 
the ray trace result and is given by 


% diff = 


computed value - ray trace 
ray trace 


The range differences were obtained by differencing the printed value of the 
range correction with the value at the previous elevation angle. Below 10* 
elevation this meant differencing range corrections over a 1 * angle interval. 

From 15° to 30° tlje differencing was over 5° angle interval, and from 
45° to 90° the interval was 15° . This exercise was not meant to simulate 
range rate corrections but rather to check for numerical stability in successive 
calculations of the range corrections. 

The heights at which 6 was computed were 1IM-G. 75YM, HM-0. 58YM, 

HM-0.5YM, HM -0. 33YM, HM, llM+150 km, H0+H1 . Ill, HUH 2 , H2 , 

2 2 

650 km, 850 km, 1500 km, and 3000 km. The quantities YM, HO, HI, and 
H2 are profile parameters as defined in figure 1. In early tests, only one 
point was used below HM and the highest point taken was at 1500 km. This 
proved to be inadequate to properly define the slopes of 6 at the end points, 
and since the spline interpolation exhibited a marked dependence on the slopes 
of 6 at these points, additional heights had to be included. Therefore the 
points at HM - 0. 75YM and at 3000 km were added into the calculation of 6 , 
for the sole purpose of allowing more accurate calculations of the first and 
second derivatives at HM-0. 58YM and at 1500 km. 

The combination of points marked 1,2, 13,9 and 10 in Tables 1,2 and 3, was 
the smallest number of points which still provided a consistently accurate 
set of range and range difference corrections at all elevations and all azimuths. 

The apparently large relative errors in the range differences seen in Tables lb, -J 

lc, and Id can be ascribed to the small absolute value of these quantities, j 

1 

i 
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in that particular azimuthal direction. In other words, percentage errors can 
sometimes be deceiving. 

The combination 1, 2, 13, 9, 10 corresponded to the seven heights HM-0.75YM, 
HM-0. 58YM, HM, HM+150, 850, 1500, and 3000. It should be noted that this 
deviation function, derived at a combination of heights which are not tied 
to the reference profile parameters, proved to be numerically more stable than 
the deviation functions whose defining heights were dependent on the reference 
profile. Furthermore it became evident in the course of this study that 
the points which most affected the accuracy of 6 were the three points 
HM-.58YM, HM, and HM+150. This indicates that great care must be taken 
to properly define the gradients along that portion of ray path which traverses 
the ionosphere within 1 150 km of the peak density. 

The optimum number of points needed for ground tracking at low elevation 
angles is seven. Since SST ray paths may be as much as 207o longer than 
such low angle cases, the optimum number of points for the longest SST 
ray paths was therefore set to nine. This seems a reasonable first guess at 
the optimum number of points to use in defining deviation functions in SST 
geometries. Further work must be done in this area to properly define the 
point function NOPF. 

B. Optimizing the Quadrature Integration 

The initial evaluations of Gaussian quadrature integration of the Bent 
ionospheric profile (fig. 1) were performed assuming a homogeneous 
ionosphere with no gradients (ref. 8) . Under these assumptions, integration 
errors of less than 0. 01% at all elevations were achieved using 3-point 
Gaussian quadrature in the bottomside bi-parabolic layer and 2 quadrature points 
each of the six topside layers (a total of 15 points). Increasing the number 
of quadrature points to 4 and 3, respectively, was sufficient to cut the errors 
to less than 0.00001%. 
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In an attempt to reduce the total number of quadrature points, the integration 
limits were modified by merging several layers. Thus, instead of integrating 
through 7 distinct profile layers, the number of integration layers was 
reduced to 3. The integration tests were again repeated using 4 and 5 layers. 

The fewer the number of layers, the more quadrature points were required 
to obtain accurate integrations. In fact, it took from 25 to 30 quadrature 
points in these tests to achieve the 0.01% accuracies obtained when integrating 
thru the seven analytic layers of the profile. Hence, it was resolved that 
all quadrature integration of the Bent ionospheric profile would be carried out 
between the limits of the ionospheric layers. As expected, in the presence of 
gradients, more quadrature points are required to maintain integration accuracy. 
The results of the numerical integrations listed in Tables 1,2, and 3 were 
all obtained using 4 quadrature points in the bottomside and 3 points in 
each of the topside layers. The range corrections are found to agree with . 
ray trace results to within 1% relative accuracy. Increasing the number of 
quadrature points to 6 and 5, respectively, had little or no effect on the 
range corrections. The changes that were noted were in no case larger than 
0.05% in the range and they provided no net improvement to the corrections. 
However, the increased number of points did add some numerical stability to 
the calculations, for the range difference error were reduced by as much as 
a factor of 2 ! The number of quadrature points used in the GTDS formulation 
of the SST refraction correction algorithm has been left at 4 and 3 pending 
evaluations using real range rate data, at which time these numbers may 
be increased as necessary. 

One can conclude from these comparisons that in the present formulation, 
it is the spline function and not the quadrature integration which sets the 
limit on the accuracy of the refraction corrections. 
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VI. CONCLUSIONS 

Comparison of the refraction correction algorithm against ray trace 
calculations confirms the hypothesis that electron density variations in the 
ionosphere can be represented as the product of a reference profile and a 
gradient function. The reference profile represents the basic vertical 
distribution of electrons while the gradient function modifies the vertical 
distribution to fit the actual density along a particular line of sight. The 
gradient function is computed, in tabular form at specified points, as the 
deviation of the electron density from the reference profile density. Up to 
twenty two such points may be used to specify the deviations. Values of the 
deviations at points between table entries are obtained by interpolation 
using a cubic spline. The reference profile for any integration path is 
defined at the point at which maximum electron density is encountered. 

Integration of the reference profile and deviation function proceeds in 
layers corresponding to the layers over which the reference profile 
maintains its analytical form. The integration is performed by Gaussian 
quadrature using 4 quadrature points in the bottomside layer of the ionosphere 
and 3 quadrature points in each of the topoide layers. 

The refraction correction algorithm has been successfully implemented into 
the GTDS, and it was used to correct SST range and range rate observations 
of the GEOS-3 satellite as relayed from ATS-6. At S-band, computed 
corrections for this orbit varied from 0.5 to 10 meter for range measurements 
and 0. 1 to 1 cm /sec for range rates. The accuracy of these values is not 
known at present, and in light of certain inconsistencies which have been 
noted in some range correction calculations, not much more can be said 
about the GEOS refraction corrections at this time. 

Ground tracking of a satellite, at low elevation angles and through regions 
of large electron density gradients was used to simulate SST geometries. 
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Comparison of algorithm calculations with ray trace results at these low 
elevations indicates that errors of less than 1% can be achieved for range 
corrections and error levels of 2% to 4% can be achieved for range 
differences. The limiting factor in achieving or improving these error levels 
is the deviation function. From ground based Tesults it is expected that 
22 points, properly distributed, should suffice to attain these accuracies 
for the worst case SST configurations. Confirmation of this will have to 
await further evalution studies. 
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Figure 5. 


A hypothetical situation showing the relationship between the profiles N and N 
and the corresponding gradient function G = N / N ref • 
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Figure 3. Definitions of some of the important points and 
parameters used in obtaining satellite-to-satellite ionospheric 
refraction corrections. Further explanations and definitions 
can be found in the glossary (Appendix H ). 
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Fig. 4.b. World Map of Height of Maximum Ionospheric Electron Density Hmax at 12 h U.T. on 
15 January 1969. Contours Every 20 km. 
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Fig. 4.c. World Map of Integrated Vertical Electron Content as Predicted by Bent Ionospheric Model 
in Units of 10 17 el/sq. meter at 12 h U.T. on 15 January 1969. 

Contours Every 5 x 10 16 el/sq.m. 
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Figure 6 . Radio propagation through the ionosphere 
R is the radio path, G is the straight line or in 


vacuo path 
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CASE 2 


igure ft. Target satellite above the ionosphere; single integration 
path from PI to P2; reference profile at PC. 
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CASE 3 


Figure 9. Target satellite in the lowest layer of 

the ionosphere (helow height of maximum 
density HM) ; single integration path from PI 
to P2; reference profile at HM. 
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CASE 4 


Figure 10 . Target satellite near occupation so that ray 

path exits and reenters the ionosphere through 
the bottom (HMYM); two integration paths (PB, PI) 
and ( PB,P2); two reference profiles at HM 
along each path. 
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Comparison of r&y trace refraction corrections vs. homogeneouR ionosphere modified by a density 
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Comparison of ray trace refraction corrections vs. homogeneous ionosphere modified by a density 
gradient function. Test case at Kashima, Japan. 12* 1 U.T. on 15 January 1969. 
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Comparison of ray trace refraction corrections vs. homogeneous ionosphere modified by a density 
gradient function. Test case at Kashima, Japan. 12 h U.T. on 15 January 1969. 
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Comparison of ray trace refraction corrections vs. homogeneous ionosphere modified by a density 
gradient function. Test case at Kashima, Japan. 12^U,T. on 15 January 1969. 
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Comparison of ray trace refraction corrections vs. homogeneous ionosphere modified by a density 
gradient function. Test case at Kashima, Japan. 12* 1 U.T. on 15 January 1969. 
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Comparison of ray trace refraction corrections vs. homogeneous ionosphere modified by a density 


gradient function. Test case at Kashima, Japan. 12 n U.T. on 15 January 1969. 
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Comparison o £ ray trace refraction corrections va. homogeneous ionosphere modified by a density 

gradient function. Test case at Kashima, Japan. 12^U. T. on 15 January 1969. 
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Comparison of ray trace refraction corrections vs. homogeneous ionosphere modified by a density 
gradient function. Test case at Kashima, Japan. 12^U.T. on 15 January 1969. 
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Comparison of ray trace refraction corrections vs. homogeneous ionosphere modified by a density 
gradient function. Test case at Kashima, Japan. 12^U.T. on 15 January 1969. 
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Flow diagram for Case Z (Figure 7} 
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Flow diagram for Case 3 (Figure 8) 
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Profile. A profile is then computed at the 
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Flow diagram for Case 4 (Figure 9) 


ORIGIN AE E5GS 19 
OF POOR QUAUCfl 





PRECEDING PAGE BLANK NOT FlUm 



LtG,S~i 














DETERMINE NOP USING* 
function nopf^M,Eca) 
WHERE W = HEIGHT OF PI 


EC A = EARTH CENTRAL ANGLE 

&u&tenwlj> &y C?i , pO 





1 

A - distance, between 


A x 0 


PO AND PI 




xU) = A 


XU) - A 


YU) * 0 


YU) « 0 


SIDE = 1 


SIDE - - l 


U - i 


N - L 

ORIGIN AD PAGE 1 

P(l) " A 


Kt) - A 

OF POOR QUAKE 

IP - i 


IP * i 


3L = 0 




JL - 0 



r<F TE R.MLNE Nop uSLNG* 

function nopf Ch jECA) 

WHERE H = \4EIGvWT Of PC 

EcA EARTH CENTRAL AN^LE 



* NOP is* USUALLY ONE LEES TM/\N HE NUMBER OF POINTS NEEDED To 
PROPERLT EVALUATE. THE ELECTRON LiiMblTT JEVTATIOM FUNCTION BETWEEN 

two Points in t he ionosphere . 










Jg M m mi M 

-m Ms,., 




E» 


SAffl 


NoP 

U&iN<« FUNCTION 
NOPF (, H , G.CA) w»it% 
H * HEIGHT OF FC 
t<-A 5 EARTH central 
A tt&U SuBTEHXD 


ssssss 


nopp - nop -OL 
loop: K=l,NoPP 

N N * 1 

Impute An* store 
Positisoms *(n) 

compute an* store 
CoRRE ^ Po M 1> 3C js» ft, 
DENSITY RATIO* Ye.**) 


SIDE - i 


*II>E. 


IP *IP ti 

PUP) = x(nop) 


fall GH.F& 

RESEQUENCES ARRAVS 
X a Nis y 


CALL "iNl 

SETS AMWVfS 
Fo(V. iptINE 
tNTtRPoUATloM 
OF DEVIATION 

function 


L 

T 


SET LINEAR 
Interpolation 
op deviation 
Function 

I 

' 

NO U«AJ>ltNTt 
PiQUIRED 


















CONVERGE, to V»» l(iMT of NAklMUM ELECT tOfsi KNitTV HM 


/is nN^ 

u^ iuov 

\'T / 


■\t OK 
HO Low 


°W< 2 ^f quS ^ 16 

0 s*>° 8 


HhPK Vim 

Cohn*Tt 1 '»i.ii,ih MAMi 
PMM OF T*l I'JU'tSi .* I loti 
OF LU»E V'ft . * t \ K[ , 
HEIGHT }«v 


Afco VI 


PRoflLf 
AT PI 


/!» 

}&o»| ML MImN 

Tat o*. 

IfcLOW 


ton en r i 

VHO; v - \ L 
t*iVi* i . t 
At i 1 1 1 *v^ 


X^Or UtiTAWCE 

tttfWKEM W-AMftPi 

Yuv : a 

N - 1 


t-OMTAiLt 
HMI?*. in ii H 

Hnj*. t 


Cowvr- #.(*(..» 


*(l) =~M "wc f | 

Kc ^ ^ tA) PC A»u> M 

YU) = *E.NMM ft At to ! 

A U-l I I 

fter *n sj it. a^pVhm i 


^ U< : l ” ‘ K ‘ nsiti,,s *« turned as follow* 
, Profile in the loop 

; ( h ;— *"»■* Uoloro.n o or Normal 
, ; ,lll ;‘; A P rw »*^‘ i* Hu*n computed at 

!; 0 , ’V M UmUuM ' lM “ r l ’’>* ThJ ratio 
U,, ‘ ! ; !>l ' ,mu>s °<1 M) ' PN(Pl) whore [>{l 
' U ‘ uj ' ttv ‘* l llu ‘ ! " 1 1 'it 1 ’t and DN{P] 

p ;V KM,S “ V l>r °‘ lu ,,Hl l,v Hu* Kelo ron 
n,!lt " ‘In- lu.ij.hi to (ho point Pl . 















LNTECjUsATOR. i 

J>fcTE<krtINt LIMITS OF INTEGRA TX.ON 



r 


**■***■■■■ -MW.: ■ grsiv. 














UtVS.M 

feTCUTTl 

AHt> VrftUQM 

F04t.TI.oH4 

POIJPOI 


HH£C0MNQ 


l*Tu*H 

1 m 


/iv rot's, 

W «lk Ma'lt 


*NT ”o 

Ktwft coo<o:nAUt or 
tKrtttt*™* 1 * TU.ru 

of\i»lt <Toi ( rot) amo 
SS»0 U»% HfelftHT 


SowuTtoH 

// 

V 


fl^ETU.O.N 


y\%, rotx, 

(M uKA^fi 


1 ORIGINAU ?At* „ 
°8 flOOR Ql/AUT^ 


n ■ tot 

P 2 -VU 


FVHfi *OF*J ®* 
lv*W«t AW40ACH 


/XI pc's. 

MTWur or >- 

■ouTXtfcE 

|(.pt,n) 

COMPUTE 

VEftUCAV. 

p^oFue 

AT PI 


TN 4 *PE 

pt.ph 


Six PiX mow 
•'Axovt 


r» VNP W«T O' ,. 

Ary^Ac't 

TO 


COMPUU 

VfcfcUC/tt. 

fKOPVCfe, 

AT rc- 


ec \ 

y l \ x. A&v'tt 

.Nl'ivvs •''*• ku^'^V ■«—■«' 

\ 14 M / 


M.VOW 


Pt = POl 
Ft * POl 


* A. Solution consists ol on*’ point ol 

intersection rtt lvoi K Ut 3500 km between 
the end points POl and FO«i. H is 
named FT V . 

H, n„M‘o are £ Into function points inside 
Hie interval (Ft>l, FO£). They »*« 
nunuut FT l and FT 3. FT l is ’ loser 
u> FOl. FT£ is closer to FO£. 

0. No intersections between FOl and F02 
hul mu’ or two intersections may exist 
outside the interval. 

Ar ovv > — >. 




; 


L£ 6 'i* 2 . 

„ 












A - DISTANCE BETWEEN 
PO And Pi. 

XU) * A 
YU) = 0 

SIDE = 1 

N - l 
Ki) 'A 
IP = i 
01 = 0 


DETERMINE. NOP USING, 

function nopfU)^ca) 


WHEB.E N s HEIGHT of PI 


EtA = EARTH central ANGLE 
£uetend&j& by (.PI , t a.) 



A x o 


XU) = A 
YU) * 0 


SIDE - -l 


N - 1 


Kl) ~ A 
IP - i 
3 L * 0 




fry 


p^ 

°e «xfe fit; 



*■ NOP is USUALLY ONE LESS than The NUMBER OF POINTS needed to 
Propfrly evaluate the electron density deviation function between 

Two Points In 1 _ He ionosphere . 

























l* hew 

»1PE 


r 

Loo 9 
I - l 

,2. ■ 



END o* 
tooP 

V i 


If* L - 1 
set pi - pv. 
If l = I 
JET PI - pi. 


/}' Pc \ 
AtaVf oft. INk u'tWS 
VhMfM s'' 


AT ‘'ft 


LEfi,S 


A * HirAMts 

*tr ftiPt 
?& = ?i 


A x 0 
VB = PC 





P 


ORIGIN AU FAGB W 
OF POOR QUAUTt 











CONVERGE TO Bottom op ionosphere HMYM 


HMYMPfc * M*YM 

PtN& COOO»M^TE6 Of 
INTEHiECTIOrt P& OF 

line (.pc >pi> amj> 
height hmym 


COMPUTE 
PROFILE 
AT pp 


y\% pl \ 

^ftOVE. fl(L *£i.Ov 


IP LOW 


AT ©<*- 
Aao\/F 



COMPAQ 

KMX HP* TO 
vmvn 


A - 

G.C tvwc. c nj 

PC \*r> PB 


couViftfeEP 
\ l X 


Y MFC. = YM 


0 rh''NVC PAGEtS' 
01 f v>oO» QOAiyfll 


F 5 



& 


fc 


g* 


H 


K$f- 

W 






*itoia 










NlER.<iE. TO height of maximum ELECTION DENSITY 


/Ss i*i\ 

AiiiWt 0^ kuo 


i\r op. 
»e u.iw 


OF Uk,t '.pft , V! and 

height Pm 


Copv^ r ( 

V E <M r ,\ u 
**«.' )i t , f 
A v v> \ i fv^ 


HMfft ,o 




ConVp Ittnfjj 

X t / 


I *U) * 1)I&TM0CE 

®* * * B *A> PC A Nil M 

YU)=&e.nmty u AtIo 

A LX) - WM AfN v r 

AN.M'WM 

YU)-0 
N :1 


tonmrt 

y eATtCAU 

•'ROFIte 



jboyj OA ft (Uhl/ 

X HU y 


At oft. 
%Eko« 


XU}= i>rwrAk»c^ 

aE rwetM pt A»*ft pi 

YU) -' o 
N - i 


I Ratio ol dons ii in * : . r 
The l ist / lor,,u ’ c ! *s follow 

. Uht °n verged profile i„ the l 

'’f’K/CNt^X^W, 

>.• rrr; pi . • ,mi w <c 









Ym-vmpb 




N - N +1 1 

kCn^A I 
0 


I 


jDETtfcttrw* Nop usiwg* 
pUwCrxoM VOPf^MjEC.^) 
\»ME<V£ W=MfIC>Hr Of Pfc 

SM- EARTH UtJTML A*nI£ 
^UBTEMDt-0 (.pa, pi'' 

J 3 JLz J L 


ORIGINAD PAT5K C 1 
OF POOR QUALITY 










INTEGRATOR £ 

J>ETfc RHINE Uj.fi ITS OF INTEGRATION 


LOOP OM 
HEIGHTS 
HTSCt) 
l - 3 ,« 


^END - 
OF Loop 

N. i v 



cow^Tt VNTta^ecrxoMfi 

PNC AND fHV 04- 
M6I<«*MT MTS Cl) 

WtTM LIME (.Pi ,? 5 .) 



I>t s DISTANCE ftETvvEEN 
PC AN© PHI 
IP = IP ♦ 1 

P(.IP)- M 


OHS 

tNftHSCUON, 


SOLUTION 


NO 

tNTC&ttCTtON 


two 

INTERSECTION* 


D 1 = MVAMCC BETWEEN 
PC and PHL 

©1 *©t&rANLE BETWEEN 
PC AN© PHI 

IP ~ IP + L 
PCtp) - -J>1 

IP ' IP T 1 

I r 











INTEGRATOR z 

PERFORM Nun EPICAL INTEGRATION OF TOTAL 
electron content /TT\ 


CALL GfcFfc 

*fc.&&Qtt£.Nfct$ 
ARRAV f> 


Loop 


I - 1 ,1? 


/ Emu x 
tooP 

\ 2 / 


RETURN 


CALL GiQUAP 

To iNTL&,ftATE 

FROM P(l^ To P(l + i) 
STotLt RESULT 
I* "PART** 

VSIT 'VNT + PART 


qriuinad PA^k 

DE POOR Qf v'.'i 








APPENDIX g 

JCL AND JOB SET-UP FOR 

SATELLITE TO SATELLITE 
IONOSPHERE RUNS 








J CL SETUP TO RUN IONOSPHERIC 
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//♦MAIf* ClAS5*OVERNITfc - 
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//SUBRLlSr 01) SVSlHJT»A» SP AC t • (CVL * 10,5I,RLSEI , ON 1 1*1 01 Sh ,91 - 
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Satellite to satellite subroutines SST and SSTCNS are modified using the 
procedure PACKVPD. Th • subroutines are members of the data set 
ORBIT. GTDS . R!>. PORT «. ordaining the source decks of the Rk D version 
of GTDS . This data sot is 1 he procedure input data set as specified by 
the DD card PACKIN. The output data set specified by DD card SEQOUT 
carries the temporary name &SOR and contains the updated routines 
SST and SSTCNS.lt is passed to the next job step where other subroutines 
are added. 


// EXEC UCL£ u ,PA*N.t 1NK=«LM, LIST, MAP, OVLY * S I 2E* ( 340K , 56K 1 • , 

// RfcChiN.t INK= 4‘>0K, KfcGI CN.GO=i»50K 

//P0SUP.0ATA5 l)U ♦ 

Other subroutines to be modified are members of the data set ORBIT. 
GTDS.SORSLIR. FORT containing the source decks of the operational 
version of GTI\5. In the first step oi procedure UCLEG they are modified 
and added to the data set &SCR. In the second step of UCLEG the 
in S»SCR are compiled and passed to the linkage editor. 
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(modifications) ° F P °° R QUALITY 

$OUTCOR 
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Distance between PC and PI. If PC lies between PI and P2 then 
A~0. 


Subroutine which performs the numerical integration of total elec- 
tron content (XNT) between points PI and P2 using Gaussian 
quadrature. GOUAl) integrates by segments having end points 
defined by the array P. 


Subroutine used to resequenee an array of up to 20 values in 
monolomcally iiu reasmg order Irom largest negative value to 
largest positive value. GRES can op. rate on a single array or 
on a pair of arrays. U one of the arrays contains positional 
values along (PI.P2), then GRFS will delete redundant points which 
are located too close ( bkm) to each other. 


Height of maximum electron density, meters above earth surlacc. 


Stores previous value of ionospheric profile parameter MM; used 
when converging on height HM. 


Height ol bottom ol ionosphere; meters above earth surface. 


Stores previous value of UMYM; used when converging on height 

HM YM . 


Number of elements in array P specifying the number of segments 
to be used to numerically integrate electron content between PI 

and P2. 


Used in integrating electron content, if PC lies between PI and P2 
then MFCS 2 and the integration segments are (PM C) and (1 G.t 
If PC lies outside (PI. 1*2) then LEGS 1 and integration segment is 

(PI.P2). 


Number of elements in arrays X and Y used in calculating electron 
density deviation function. 


Number of segments dividing the interval (PI. Pi) into equal earth 
central angles. Used in determining the number ol points N to 
be used in defining the deviation function which describes the 
variations in electron density between PI and P2 or any othei 
integration path. 
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P 


PB 


PC 

PI 


(j) P01, P02 


PT 1 , PT2 

P1,P2 

SIDE 


SPLN2 



Function which upon entry returns the number of points needed 
to properly define the electron density gradients between any two 
points in the ionosphere. The input parameters are a) the earth 
central angle subtended by the Line segment joining the two points 
and b) the height above the surface of the earth of the point of 
closest approach of the Line segment. 

Array of positions along (P1,P2) marking the end points of the 
segments used to numerically integrate the total electron content 
between PI and P2. The points chosen are the intersections of 
(P1,P2) with the heights HMYM, IlM and the five heights 
delineating the topside exponential layers. 

Coordinates of intersection of the line of sight and the bottom of 
the ionosphere. 

Point of closest approach of the line (P01.P02) to earth center. 

Takes on value of coordinates of points PI or P2, as appropriate. 
Used to define line segment end point for integration, either 
(PC, PI) or (PC, P2). 

Body fixed XYZ coordinates of end points of line of sight. 

Ground to satellite mode: 

P01 is the position vector of ground station 
P02 is the position vector of satellite. 

Satellite to satellite mode: 

P01 is the position vector of target satellite 
P02 is the position vector of relay satellite 

Body fixed XYZ coordinates of intersections of line of sight 
and 3500 km sphere. 

Body fixed XYZ coordinates of end points of that segement of the 
line of sight (P01, P02) which lies entirely within the ionosphere. 

Parameter used to distinguish between path segments on opposite 
sides of PC. If PC lies between PI and P2 tlxen for segment 
(Pi, PC) SIDE= - 1 and for (PC,P2) SIDE=+1. If PC lies outside 
(PI, P2) then SIDE =+l. 

Subroutine which computes and sets up arrays necessary to perform 
a spline interpolation within an array of vlues. Interpolation is 
performed by subroutine SPEIN. 

Array of positions along (Pl,P2) at which the density ratios Y are 
computed. ^The positxofxs~a T'e" irT meter?* 61 dfstan'Ce from PC 
along (P1,P2). 
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XNT 


Total electron content integrated along the line of sight from 
P01 to P02. 


Y 


YMPB 


Array of electron density ratios used in computing the deviation 
function which defines the electron density gradients along (PI PZ) 
(PI, PC) or (PC, P2) depending on the geometry, 

DEN (A. , 9 , h) 

DENREF(h) 

Where DEN indicates electron density at position \ , & , h and 
DENREF is density of reference profile at height h. 

Ionospheric profile thickness parameter YM of lower ionosphere 
at the point PB. 
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